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Meso-scale models and approximate solutions for solids 

containing clouds of voids 

V.G. Maz'ya* A.B. Movchan^ and MJ. Nieves^ 


Abstract 

For highly perforated domains the paper addresses a novel approach to study mixed boundary value 
problems for the equations of linear elasticity in the framework of meso-scale approximations. There are 
no assumptions of periodicity involved in the description of the geometry of the domain. The size of the 
perforations is small compared to the minimal separation between neighbouring defects and here we discuss 
a class of problems in perforated domains, which are not covered by the homogenisation approximations. 
The meso-scale approximations presented here are uniform. Explicit asymptotic formulae are supplied with 
the remainder estimates. 


1 Introduction 

Meso-scale approximations have been introduced and rigorously studied in iTlSl l22l l24l . Physical applica¬ 
tions in composite systems in electromagnetism were also addressed in the earlier papers |IZl|8l. The study 
of Green's kernels as well as asymptotic analysis of solutions to eigenvalue problems for dense arrays of 
spherical obstacles was performed in 1281 . Compared to classical homogenisation approaches (see l3l l29llT0l ), 
the meso-scale approximation does not require any constraints on periodicity of the microstructure, and it is 
uniformly valid across the whole domain, including neighbourhoods of singularly perturbed boundaries. 

Prior to the development of the meso-scale asymptotic approach, many papers and monographs (see, for 
example, ||5j |6l HlJ IT2l l have appeared which model singular perturbations of various domains. Examples 
include domains with irregular boundaries, thin components or domains containing either a single small 
defect or several defects. The method of compound asymptotic expansions of solutions to such problems is 
described in l26l |27| . In particular, for domains with small defects, asymptotic approximations have proven 
to be superior to the finite element method (FEM), even when the overall number of defects is chosen to 
be large 1191 . For domains with perforations, the approximations presented in Il26ll27ll use model problems 
posed in the domain without defects and problems posed in unbounded domains, in the exterior of individual 
inclusions. Integral characteristics of the defects are used here in cormection with the energy of model fields 
in the exterior domains. For rigid inclusions we refer to the capacity of the inclusions, whereas for voids we 
use the dipole matrix, that correspond to the Dirichlet and Neumarm boundary conditions, respectively. 

The method of compound asymptotic expansions has also led to the development of uniform approxi¬ 
mations for Green's kernels for domains with small defects for the Laplacian, corresponding to a variety of 
boundary value problems involving rigid inclusions IITtIITsI , voids IITTII and soft inclusions 12111 . Approxima¬ 
tions for Green's kernels in long rods have also appeared in 11611 . There also exist several approximations for 
Green's tensors of vector elasticity for solids with rigid inclusions ll9l l20l and holes with traction free bound¬ 
aries 1231 . Meso-scale approximations of Green's function for the Laplacian in a solid with rigid boundaries 
has been derived in HSl . 

A systematic presentation of the theory of meso-scale approximations in densely perforated domains is 
given in the recent monograph l23l . In particular, it was demonstrated that uniform meso-scale asymp¬ 
totic approximations are of high importance for the analysis of fields in solids containing non-uniformly 
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The collection of voids 1 < i < ^ 


Figure 1: The solid Dpj containing a cloud cv of voids. 


distributed clouds of small voids or inclusions. In such configurations, the traditional computational ap¬ 
proaches like FEM are inefficient. 

Recently, the method used to develop meso-scale approximations for scalar problems posed in solids with 
many small voids and inclusions has been extended to the Dirichlet problem of elasticity in solids with a 
cloud of rigid inclusions Il24t . In the present paper the approach of 1231 is applied to a mixed boundary value 
problem of vector elasticity in an elastic solid, which contains a cloud of many voids whose boundaries are 
traction free. The number of voids is denoted by N ^ 1. Each void is a concentrator of stress, and analysis 
of boundary layers is carried out in terms of special classes of dipole fields, which characterise the shape of 
voids and elastic properties of the material. The schematic representation of the porous solids with a cloud 
of N voids is shown in Eig. Two small parameters are introduced as the normalised diameter of a void and 
the minimal distance between neighbouring voids within the cloud. 

Let n be a bounded domain in IR^ representing an elastic solid. Contained in O are many small voids, 
1 < i < N whose diameters are characterised by the small parameter e and that occupy a set a; C O 
representing a cloud of voids. The sets O and j = 1,... ,N, are assumed to have smooth boundaries. In 
addition, the minimum distance between the centres O , 1 < A: < N, of each void is cormected with another 
small parameter d. The geometry of the elastic solid with many small perforations will be described by the 

set On = n\ ujlj 

In the framework of vector elasticity, the Lame operator and the operator cormected with the application 
of external tractions will be denoted by L(Vx) and T(Vx), respectively. 

The displacement field un satisfies the governing equations of static elasticity: 

h(Vx)u]v(x) -F f(x) = O , x e , 
un(x) = O , x e 30 , 

T„(Vx)un(x) = O , x e , 1 < ; < N . 
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( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 







In 11.11, f G Lco{0,n) is a vector function representing the action of body forces inside fhe perfurbed solid, 
fhaf safisfies fhe consfrainfs co n supp f = 0 and disf(supp f, do;) > C, wifh C being a posifive consfanf 
independenf of e and d. 

The consfrucfion of fhe approximafion for presenfed here depends on several model fields: 

1. fhe solufion u of fhe problem in O wifhouf any voids, 

2. fhe regular parf H of Green's fensor in O, 

3. a mafrix funcfion fhaf solves a Neumarm problem in fhe exferior of fhe scaled void whose 

(k) 

columns are known as fhe dipole fields for fhe elastic void; a re-scaling is applied fo obfain ' for fhe 


small void co 


[k) 


4. a consfanf mafrix called fhe dipole mafrix of fhe scaled void fhaf characferises fhe void's 

shape and fhe elastic maferial properties. The dipole mafrix for fhe small void cug is consfrucfed 

from by re-scaling. The geomefry of fhe voids is assumed fo be chosen so fhaf fhe maximum and 

(k) (k) (k) 

minimum eigenvalues A^ax and respectively, of fhe mafrix — ' safisfy fhe inequalities 


Cl < A**"^ 

— min 


and 


^max S ^2 £ , 


(1.4) 


for k = 1,... ,N, where Ci and C 2 represenf differenf posifive consfanfs. 


For convenience of nofafions, we also use fhe vector E of normalised elasfic sfrain, corresponding fo fhe 
displacemenf field u, so fhaf E(u) = E(Vx)u, where E is fhe linear mafrix differenfial operafor. 

The consfanf vector V and mafrices M and S are also used in fhe approximafion for ujv: 


V=((E(V,)^u(x))^ 


(= 0 ( 1 ) 


.(S(V,)^u(x))^ 


(=o(to 




M = diag{Mp^...,M^^^} 


and 


S = 


E(V,)^(E(Vy)^G(y,x))^ 


x=0(') 

y=0(l) 


if i 7 ^ i > 

ofherwise, 


where ©6x6 is the 6x6 null mafrix; also in fhe fexf below I^xn will sfand for fhe n x n idenfify mafrix. 

The main result of fhis arficle is fhe uniform asympfotic approximafion of fhe displacemenf field ujv, as 
presenfed in fhe following fheorem: 


Theorem 1 Let the small parameters e and d satisfy the inequality 

e < cd, 

with c being a sufficiently small constant. Then the approximation for u^ is given by 

N 

un(x) = u(x) + ^ |q®(x) - (E(Vz)^H(z,x))^M® 


! = 0(*) 


}cW+Rn(x) 


where C = ((C^^i)^,..., solves the linear algebraic system 

- V = (l6Nx6N + SM)C , 
and for the remainder R^, the energy estimate holds 

tr{cr{Risi)e{Risi))dx < Consf ||E(u) 

Here Consf in the above right-hand side is independent of e and d. 


(1.5) 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 
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This representation | |1.6| is uniform and it engages several classes of model fields, which are independenf 
of fhe small paramefers e and d (also see 112^ 1. 

The sfrucfure of fhe paper is as follows. Main nofafions are infroduced in Secfion|^ Model problems used 
fo approximafe ujv are infroduced in Secfion|^ The formal approximafion of u^ is fhen provided in Secfion 
1^ This approximafion relies on fhe solufion of fhe algebraic sysfem (1.7l and fhe solvabilify of fhis sysfem is 
sfudied under fhe consfrainf in Secfion 1^ Then, in Secfion fhe energy esfimafe < |1.8| for fhe remainder 
of fhe approximafion is proved. Simplified asympfofic approximafions for U]v are fhen given in Secfion 
Following fhis, conclusions and discussion are given in Secfion Appendix A confains a local regularify 
esfimafe used in fhe proof of fhe energy esfimafe < |1.8| . In Appendix B, a defailed proof of infermediafe 
sfeps used fo show fhe solvabilify of •HZI is presenfed. Finally, in Appendix C, we show fhaf for cerfain 
geomefries, dipole characferisfics can be consfrucfed in fhe closed form for fhe case of spherical cavities and 
explicif represenfafions are given. 


2 Geometry of the perforated domain and main notations 

A domain O c will be used fo denofe fhe sef corresponding fo an elastic solid wifhouf holes, wifh smoofh 
frontier 30. For a small positive paramefer e > 0, fhe open sef is defined in such a way fhaf if confains an 
inferior poinf has smoofh boundary and a diamefer characferised by e. The collecfion of sefs 
1 < 7 < will represenf fhe small voids confained inside fhe sef O fhaf are subjecf fo some furfher geomefric 

consfrainfs discussed below. In fhis way, we define fhe perfurbed geomefry On = n\ If is also 

assumed fhaf a small paramefer d characferises fhe minimum disfance befween poinfs in fhe array 
and fhaf fhis minimum disfance is 2d. Anofher geomefric consfrainf is fhe assumpfion of fhe exisfence of a 
sef CO fhaf safisfies 


N 


N 


C CO, disf( y co^'^ 
j=l /=! 

If is also useful fo infroduce fhe mafrix functions: 


S(x) = 

and 

?(x) = 

These mafrices satisfy fhe condifions 




,dco) > 

2d 

and 

disf(3a;, 30) > 1 . 

2-1/2 

X2 




2-1/2 

Xi 

0 


1 (2.1) 

! 0 



2-^/^X2 ) 



2- 


0 ^ 

1 



0 

2-1/2x3 

(2.2) 

0 

-2 


-2-^/^X2) 

1 


S(Vx)^S(x)=l6 


x6 


S(Vx)^?(x) = 06 


x6 


where fnxn and O^xn are fhe n x n idenfify and null mafrices, respecfively. 

The mafrices ^ and E also lead fo a compacf form of fhe firsf-order Taylor approximafion for a vecfor 
funcfion u abouf x = O 


u(0) = ?(x)?(Vx)^u(0) +E(x)S(Vx)^u(0) 
and allow fhe Lame operator L(Vx) to be defined as 

L(Vx) :=E(Vx)AS(Vx)^, 

wifh 


■o(b 


( B ©3x3 \ 

(\ + 2}i 

A 

A \ 

A = 

B = A 

A + 2ji 

A 

VOsxS 2^(13 / 

V A 

A 

\ + 2^J 
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The corresponding traction operator Tn^Vx) is then 

T„(Vx) :=E(n)AE(Vx)^, 

which will be applied on the boundary of an open set with n being the unit-outward normal to the set. 

The strain tensor e(v) = stress tensor cr{v) = [o'ij{v )]^and the tensor of rofations i]{v) = 

['?//■ ^ vector field v take the forms 


and 


e(v) = 2 ((^ ®v) -I- (V (g) v)'^) , £r(v) = Atr(e(v))]t 3 -|-2pe(v) , 


»/(v) = -((V®v)-(V®v)^). 


The mafrix J = where /i'i is fhe column of fhis mafrix, is 

/ 0 -X 3 X 2 

J(x) = X 3 0 -xi 
\ —X2 X-[ 0 

and fhis plays a role in fhe descripfion of fhe overall moment acting on an elastic body. It is noted that 


(2.3) 


'0 0 0 ' 


0 0 -1 


0 1 o' 


^(/(i)) = 0 0 1 , »/(/(2)) = 0 0 0 and =-100 


0-10 


1 0 0 


0 0 0 


The strain and stress vectors denoted by E and N, respectively, are defined by 

E = (fill, £22, £33, V^£12, V^£13, V2£23y , 

N = ((711,(722, ( 733 ,-\/2cr32, •\/2cri3,-\/2(723)^ , 
and can also be infroduced fhrough fhe mafrix operafor •ED as 

E(v) = E(Vx)v N(v) = AE(Vx)v, 

for a vector function v. Nofe also fhaf fhe quantify S(U) = fr(e(U)e(U)) can also be represented as 

I3 03 


S(U) = E(U)^ 


,03 2-II3 


E(U). 


(2.4) 


(2.5) 


( 2 . 6 ) 


3 Model fields 


In fhis secfion, we discuss fhe model fields used in fhe meso-scale approximafion of U]v in defail. We begin 
wifh an infroducfion of fields defined in fhe unperfurbed sef O: 

1. The solution of the exterior Dirichlet problem. The vector field u is a solufion of 

L(Vx)u(x)-hf(x) = O , xeO, (3.1) 


u(x) = O , X e 30 , (3.2) 

where f satisfies fhe same condifions as in fhe sfafemenf of problem (1.1 1-( |1.3| , and fhe same nofafion 
will be used to represenf fhe exfension of f by zero inside fhe voids (X’Fh 1 < j < N. 
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2. The Green's tensor for the solid O. The notation G will refer to the Green's tensor in the domain O, that is 
a solution of 

L(Vx)G(x,y)+^(x-y)l 3 = OI 3 , x,yen, (3.3) 

and safisfies fhe homogeneous Dirichlef condition 


G(x,y) = OI 3 , X e dCi,y e O . (3.4) 

The regular part H of fhis fensor is represenfed by 

=T(\,j) - G{\,y) , (3.5) 


where T is fhe Kelvin-Somigliana fensor 


r(x,y) 


1 1 

Sn}i{\ + 2}i) |x — y| 


(A + 3/r)]t3 + (A + }i) - 


y) 0 (x - y) ] 

|x-yP j' 


(3.6) 


and 

L(Vx)r(x,y)+<5(x-y)l3=0l3. 

The above problem then implies that H(x,y) = (H(y, x))^, x,y e O. 

Next, we introduce the boundary layer fields for fhe small voids, known as fhe dipole fields IIT7ll25l : 


3. The dipole fields for the voids. In fhe consfrucfion of fhe boundary layers for fhe asympfofic algorifhm in 
the vicinity of fhe void fhe physical fields known as dipole fields, will play an essenfial role. They 
are defined as funcfions of fhe scaled variable = e“^(x — O^i^) oufside of fhe scaled sef = {^j : 

G The dipole fields form the columns of fhe 3x6 mafrix where 


L(V^JqW(?^.) = C)3x6, e ]R^a;0■) , (3.7) 

T„(V^JqW(?^.) = S(nW)A , e amO') , (3.8) 

qW(^^.)^ 03 x 6 , as (3.9) 

where n^t) is the unif outward normal fo 

The righf-hand sides in fhe Neumarm boundary condition < |3.8[ are subjected to the constraints that the 
total force on boundary and fhe resulfanf momenfs are zero: 


Lu) = 03x6 , (3.10) 

Lu) = 03X6 . (3.11) 

A special mafrix mOI, wifh consfanf enfries, is also required fo consfrucf fhe leading order behaviour of fhe 
mafrix qO) af infinify and fhis is called fhe dipole mafrix. The behaviour of far away from fhe void 
is described in fhe nexf lemma (see nzma). 

Lemma 1 For \^j\ >2 the matrix qOI admits the form 

= -(s(v^yr(^^,o))^MW +o(|^^.|-') ■ 
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4 Formal meso-scale approximation for ujv 

In this section, the derivation of the meso-scale asymptotic approximation for ujv in Theorem is formally 
derived. Firsf we nofe fhaf in what follows, we will need fhe mafrices (x) = i^j) and = eSmO'). 
According to IIT7ll25] , the dipole matrix is symmetric negative definite. 

In the next Lemma and the following texf fhe nofafion Consf will represenf differenf posifive consfanfs 
independenf of fhe paramefers e, d and N. 

The meso-scale approximafion for fhe displacemenf field ujv is now defined by 

Lemma 2 The formal approximation o/ujv is given in the form 


N 

un(x) = u(x) + (Qe (x) - (S(Vz)^H(z,x)) 
k=i 




! = 0('=) 


}cW+Rn(x) (4.1) 


where the coefficients satisfy 


E(Vx)^u(x) ^ S(Vx)^(E(V.)^G(z,x))^M, 




l<k<N 




z 

x=oW 


(4.2) 


for 1 < j < N. The remainder is a solution of the boundary value problem for the homogeneous Lame equation in 
Cif^, with the mixed boundary conditions: 

Rn(x) = (p{\) on 30; and Tf,(Vx)RN(x) = onx E dcoP ,1 < ; < N, 


where the right-hand sides satisfy the estimates 


and 


|^(x)| < Consf 


ti|x-OW|3' 


x E 30 , 


|^(^'^(x)| < Consf 


^e(l -h e^|cW|) 


E 

k¥^i 

\<k<N 


e4|cW| \ 

x-OW|4 ) ' 


X E 3a;P , 1 < 7 < N , 


(4.3) 


(4.4) 


and (p^i\ 1 < 7 < N, fulfil the orthogonality conditions 


Jda}f 


Sx = o , [ J(x - oW )</>Wdsx = o , 1 < 7 < N . 


/dco. 


(4.5) 


Proof. The orfhogonalify condifions 14.51 follow from | |4.1| , fhe Beffi formula and fhe model problems 
infroduced in Secfion 2. 

According fo problem 1, Secfion 2, the vector function 


= un(x) -u(x) , 


(4.6) 


satisfies fhe homogeneous Lame equation for x G O^. Since bofh un(x) and u(x) safisfy fhe homogeneous 
Dirichlef condifion on 30, fhen R^^(x) = O, for x E 30. Nexf consider fhe fracfions of fhe Rj.]^ on d(jo^\ 
This condition, using Taylor's expansion about x = O^t) takes the form 


T„{V^)r\^'’ = T„(Vx)(un(x) - u(x)) = -T„(Vx)u(x) , 


= -T„(Vx)u(x) 


x=oW 


-h 0(e) , x G dcvg'^ , 1 < 7 < N . 


(4.7) 
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An approximation for is then sought as 


N 

Rn W = +Rn(x) . (4.8) 


The goal is now to determine the vector coefficients 1 < A: < N, to complete the formal approximafion. 
If is nofed fhaf fhe remainder in ( |4.8) is a solufion of 

L(Vx)Rn(x) = O , xeOjv, 


and from fhe boundary condifion for regular parf H of Green's fensor (see <3.3 1-< |3(5) ), fhe exferior Dirichlef 
condition for R]v/ is 


Rn(x) = (x)-(E(V.)^r(z,x))W)[^^^Jc 




N ,4 


= o ( r " 


, X G 30 , 


(4.9) 


where Lemma [^has also been used. In order fo derive fhe vector coefficienfs C^i\ 1 < / < N, fhe fractions 


on fhe inferior boundaries for Rjv should be considered. For x G dco£^\ according fo <4.7 
T„(Vx)Rn(x) = -T„(Vx)u(x) - T„(Vx)qF\x)cW 


- E T„(Vx){qF(x)-(E(V.)^H(z,x))^M, 
l<k<N 




!=oW 


jcW 


+0(e) +0(e3|cW|) , X G A < j < N . 

Condifion <3.8' and Lemma [^fhen provide a simplified form of fhe above fraction condifion on dcojJ'^: 


,(Vx)Rn(x) = -S(nW)A|s(Vx)^u(x) 


X=o0') 

■ E(Vx)^(E(Vz)^G(z,x))^M® 

l<k<N 


£ = OW 


cWj 


+ 0 (£) + 0 (e 3 |cW|)+ 0 ( 


l<k<N 


e4|C«| 

lx - OW| 


, X G dcdP , 1 < 7 < N . 


Applying fhe Taylor expansion once more abouf x = gives 
'f’n(Vx)RN(x) 

= -S(nW)A{s(Vx)^u(x)[^^^.^+cW+ £ E(Vx)^(S(V.)^G(z,x))^M« 


l<k<N 


2=ow'- / 

X=o0') 


+0(£)+0(e3|cW|) + o( ^ 


i|C«| 


|O0')-OW|4^ 

l<k<N 


, X G 3mp^ , 1 < 7 < N 


(4.10) 


Thus, we can remove fhe leading order discrepancy in fhe preceding boundary condifion by allowing fo 
safisfy fhe sysfem of equations 

E(Vx)^u(x) +CW + E E(Vx)^(S(V.)^G(z,x))^mW .= 0 ( 9 ^^'^ = O , (4.11) 


l<k<N 


=oW' 

x=o9) 


for 1 < 7 < N. Combining <4.6 1 , <4.8 1 and <4.9 H 4.111 complefes fhe proof of fhe lemma. 


□ 
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5 Algebraic system for and its solvability 


Before presenting the energy estimate for the remainder R]v, the solvability of the algebraic system (4.21 is 
discussed in this section under the constraint that e < cd. We first introduce some notations to simplify the 
analysis. 

Using the following vectors. 


C=((cW)^.(CW)^)^ and V = ((E(Vx)^u(x))^ .(E(Vx)^u(x))^ 




and the 6N x 6N symmetric matrices: 


M = diag{Mp^...,My'''’} , 




s = 


E(Vx)^(S(Vy)^G(y,x))^ 


©6x6 / 

the equations (|4.2| can be written as 


x=o(^) 

y=0('=) 


if i ^k, 
otherwise . 


- V = (l6Nx6N + SM)C . 


5.1 Solvability of the algebraic system jS.lf 

Here, a result concerning the solvability of the system is proved: 
Lemma 3 Let the parameters e and d satisfy the inequality 

£ < cd, 


where c is a sufficiently small constant. Then, the linear algebraic system (5.1) is solvable and 

N N 

£ |CW|2 < Const |E(u(x))|2 
M 7=1 


c=o0') ' 


where the strain vector E(u(x)) is defined in (2.4). 


Proof. By taking the scalar product of 0 with MC and using the Cauchy inequality we deduce 
(-MC,C) - (MC,SMC) = (MC,V) < (-MC,C)^/^(-MV, . 

Note that the term (MC, SMC) admits the form 


N 


MQSMC) = ^(mPcW)^- S(Vx)^(E(Vy)^G(y,x))^ . 


7=1 


^^7 

l<k<N 


x=o0') ' 
y=0(*') 




In Appendix B, it is shown that ( |5.5[ satisfies 

|(MC,SMC)| < Const d"3(MC,MC) . 
Returning to ( |5.4[ , this can then be used to establish that 

\i/2 > /_Mr r\i/2 _ ( MC,SM C) 


-MV,V 


> (-MC,C 

> (-MCC)^/^ - Const d 


We note that 


-MC,C)1/2 

_3 (MC,MC) 

(-MC,C)1/2 ■ 

(MC,MC) = (-MC,-MC) < Const max A®x (-MC,C) 

l<i:<N 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


(5.5) 


(5.6) 


and since the eigenvalues of the dipole matrices — Mg*^^, 1 < k < N are O(e^) according to (1.41, it follows 

(-MV,V)1/2 > (1-Const e3d-3)(-MC,C)i/2 . 

Estimate (|5.3| now follows from (|5.2[ and (|2.5[. The proof is complete. □ 
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6 Energy estimate for the remainder Rjv 


With the formal meso-scale asymptotic approximation of un in place, fhe energy esfimafe for fhe remainder 
ferm R]^ in Theorem is now obfained. 

Lemma 4 Let the parameters e and d satisfy the inequality 

e < cd 

where c is a sufficiently small constant. Then the remainder term satisfies the energy estimate 

j fr((r(RN)e(RN))dx < Consf +e®d“^|||E(u)|||^(j^j , (6.1) 

where the constant in the right-hand side is independent of e and d. 

Prior fo fhe proof of Lemma and Theorem 1 we infroduce several auxiliary nofafions. 

6.1 Auxiliary functions 

In fhis parf of fhe proof, auxiliary functions will be infroduced fhaf will allow fhe remainder Rj^/ fo be 
esfimafed. The firsf funcfions fo be considered are fhe cuf-off funcfions whose supporfs are locafed in fhe 
vicinify of fhe boundaries of 0 ^. 

Namely, fhe cuf-off funcfion Xe^'^ G ^ ^ k < N will be used fhaf is equal fo 1 inside fhe ball ■ 

A cuf-off funcfion Xo is also required and will allow for cerfain domains of infegrafion fo be concenfrafed 
near fhe boundary 30. Wifh fhe sef V^ = {xGn:0< disf(x, 30) < <5} we define xo G where 

V = Vi/ 2 . The funcfion xo is equal fo 1 on V 1 / 4 , and zero when x E 0\V. 

Now vecfor funcfions A: = 0,1,..., N, are infroduced fhaf safisfy fhe conditions: 


Yo(x) = -Rn(x) , 

for x E 30 , 

(6.2) 

and 

for x e dco^:^^ , p = 1,... ,N . 


T„(Vx)Yp(x) = -T„(Vx)Rn(x) , 

(6.3) 

Such funcfions will fake fhe represenfafions 



Yo(x) = E{qP(x)-(E(Vw) 

^r(w,x))^Mp') jcW 

(6.4) 


and for 1 < k < N 

Y,(x) =u(x)-^(x-OW)(?(Vx)^u(x)) 
-(D(Vw)’'H(w,x))^M®ci*') 


x= 0 (^ 
w= 0 (*^) 


-E(x-0«)(E(Vx)^u(x)) 


x=oW 


+ E {qF^(x)-(S(Vw)^H(w,x))^mF^| (JcW 


1</<N 


- E E(x-0«)E(Vx)^(E(Vw)^G(w,x))^Mp')cW 


\<j<N 


x=oW ■ 

w=Oij) 


(6.5) 


Wifh fhese choices for fhe funcfions Yj^, 0 < A: < N, if can be verified fhaf fhey indeed safisfy (6.2 1 and (6.31. 
Also nofe fhaf for A = 1,..., N if can be checked fhaf 


f (k) In (Vx)Yj;(x)dSx — O 

Jdojf ’ 


and 


Jdw, 


,^J(x-0W)T„(Vx)Yfc(x)dx = 0. 


( 6 . 6 ) 
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In the sequel, we also use the same notation Rj^/ to denote the extension of the remainder into the regions 
1 < A: < N, similar to 1301 . 

Later, the constant vectors 


rW = 



l<k<N, 


(6.7) 


and 

= ^ / (RN(x)+J(x-OW)r«)dx , l<k<N, (6.8) 

will also be required. Using these constants, a rigid body displacement can be constructed in the form 
Rn(x) + J(x)r(^) fhaf safisfies 


and 


/ »/(Rn(x)+J( 

'^^3£ 

/ (Rn(x)+J(x- 


x - 0(''))r(''))dx = ©3x3 , 

(6.9) 

0 W)rW = 0 . 

(6.10) 


6.2 Estimate for the energy in terms of the functions 

Here if is shown fhaf 


N 




fr(c7-(R]v)e(RN))dx < Consf < / |Yo|^dx+ / |E(Yo)|^dx+^ \E{'f,,)\^d\'> . (6.11) 


k=l-'°3e 


Firsf, sef 


N 


W = Rn + XoYo and V = RN+Y^X^P^k- 

k=l 


( 6 . 12 ) 


Nofe fhat according fo (6.2 1 and (6.3 , W = O for x e 30 and Tn(Vx)U = O, for x e As a resulf, 

affer applying Beffi's formula, it is possible to show that 


/ n,\ 


tr(£r(W)e(U))dx = — / W-L(Vx)Udx. 


/ Clf. 




Recall the supports of fhe cuf-off funcfions xo and Xe A k = 1,... ,N do nof infersecf, and Rjv safisfies fhe 
homogeneous Lame equafion in Cij^. Thus affer replacing U and W with their definitions in ( 6.12| , the 
preceding identity reduces to 

/ tr{cr{R^ + xo^oMRN+Exi’'^^k))dx=-'£ / ^RN-L(Vx)(xfY;t)rfx (6.13) 

which can be further simplified by expanding fhe leff-hand side using fhe linearify of fhe sfress and sfrain 
fensors fo give fhe inequalify 


/ Ox 


fr(cr(RM)e(RN))dx < Ui + £2 + U3 


(6.14) 


where 


Ui = 


E 2 = 


u, = 


tr{p-{xo'¥o)e(RN))dx 

t L,-Xk^^N-L{V.){xP^k)dx 

V —1 'A D'X^ \(^F 


k=l •"’3E \“'£ 
N 


E / w,-wtr(cr(RN)e(,y®Yfc))dx 


(6.15) 


Nexf, fo derive (6.11 1 , Uy, j = 1,2,3 are esfimated. 
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6.2.1 Estimate for Ei 

The term Zi, by the Cauchy inequality and the Schwarz inequality, admits the estimate 

< (^^tr(cr(xoYo)cr(xoYo))dx^ (^J^S{RM)d> 

Here, the quantity S(U) is defined in | |2.6[ . Since the inequalities 

Cl 1 n / (3A + 2^)2 if 0 <v<l /2 

tr(£r(v)cr(v)) < CiS(v) where ci = | ^^2 if - 1 < v < 0 


(6.16) 


(6.17) 


and 


-K / if0 <v<l/2 , 

trWvjeW) > C 2 S(v) where if - 1 < e < 0 , 


hold for a vector function v, it can then be asserted from <6.17}, <|6.18 1 and <6.16} that 


Si < Const 


S{xo^o)dx 


1/2 


’Cif, 


tr(£r(RN)e(RN))dx 


1/2 


(6.18) 


(6.19) 


6.2.2 Estimate for E 2 

Note that 


/ ^(J(x)r« - Rw^*^)) ■ L(Vx)(^®Y,)dx = 0 , 


( 6 . 20 ) 


where the definitions of r(*^) and Rj.^ are found in <6.7 1 and < 6 . 81 . Identity <6.20 1 appears as a result of the 
application of the Betti formula in \cOe^'^ as follows: 

■ L(Vx)(^F'Ffc)dx = - ■ L(Vx)(J(x - 0«)r« - R^^^^Vx 


+ 


L 


{(J(x - 0«)r« - Rn^'^^) ■ T„{V,){x^’^^^k) - ■ T„(Vx)(J(x - 0«)r« - RN^"^)}dSx . 




a(B(9\a;<9) 






( 6 . 21 ) 


The first integral on the right is zero since all rigid body displacements are solutions of the homogeneous 

(k) 

Lame system. They also produce zero traction and this together with the definition of Xe k < N, shows 

that 

“ 0 «)r« - Rn^*^’) ■ L(Vx)(;tf'r/c)rfx = /g^w(J(x - 0 «)r(*^) -Rn^'^^) ■ T„(Vx)Y,dSx , 

and owing to < | 6 . 6 } the right-hand side is zero. 

In addition to < 6.20} , the next identity is also true 

/ ^(Rn(x) + J(x - 0«)r« - R^^'^’) ■ L(Vx)(Tf (J(x - OW)^.® - %))dx = 0 , (6.22) 


where similar to <6.71 and <16.8} 


= ^ /wJ(Vx)Y,(x)dsx, l<k<N, 


and 


71= ' 


d(^) 

^3£ 


■ / (Y;t(x)+J(x-0«)l/>«)dx , l<k<N. 

JbX 
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Here (6.221 follows from applying fhe Beffi formula inside and making use of fhe facf fhaf Rjv is a 


solufion of fhe homogeneous Lame equafion in Ojv and fhaf if safisfies fhe condifions (4.51. 


Therefore, fhe ferm E 2 in ([6.15}, in combinafion wifh (6.20| and (6.221, is also written as 


2:2 = 


N „ 
k=l ■' ^3£ 


(Rn(x) + J(x - oW)r« - R«) ■ L(Vx)(xf (Y,(x) + J(x - 0«)i/.« - Y,))d> 


(6.23) 


The Schwarz inequalify followed by fhe Cauchy inequalify shows fhaf E 2 is majorised by 

Consf /m I(Rn(x) + J(x - - RN^*'^)|2dx^ 

VJ :=1 / 

/ N \ 1/2 

X E /„)mVx)(^f(Y,(x)+J(x- 0 «)i/.«-Y^))| 2 dx , 

\J :=1 / 

where R^ has been smoofhly exfended inside Then Poincare's inequalify shows fhaf in ^ = 

1.N, 


|(Rn(x) + J(x - 0«)r« - R^^)\^dx) ' < Consf 6 ( |V(Rn(x) + J(x - 0«)r«)|2dx') 


\ 1/2 


(6.24) 


Nexf as a resulf of condifion (6.91, fhe Friedrichs inequalify can be used, similar fo ||9l/ to give the estimate 


|V(Rn(x) +J(x - o(''))r(*'))|2dx^ < Const e S(RN)di 

This argument together with (6.181 and (6.23| shows that 


1/2 


(6.25) 


N 


1/2 


E2 < Const e E \ ik) ti'(‘^(RN)e(RN))dx 

\k=i / 

X (e /(„|B(Vx)(xf (Y,(x)+J(x-0«)t/>«-Y^))|2dx) 

Xk^i-'hs ) 


1/2 


(6.26) 


By computing derivatives and taking into account the definition of the cut-off functions Xki k = 1,..., N, an 
estimate for the second integrand on the right can be established in the form 

|B(Vx)(^F(Y,(x) +J(x)t/;W - Y,))|2 < Const | V(Yfc(x) + J(x - 0«),/-«)|2 

+e-2|Yt(x) + J(x - oW)i/>W - Y^|2} , (6.27) 


where the L(Vx)Y;t = O/ for x e B® has been used. 

Thus, ( 6.27| together with the application of the Poincare inequality and the Friedrichs inequality inside 
Bg^^ leads to 


/ N \ 

E /(„|B(Vx)(xf (Y,(x)+J(x-OW)t/.W-Y,))|2dx 

< Const e"i ( ^ [ S(Yfc)dx 
\k=i 


(6.28) 
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Combined with ( 6.26} and the fact that 

(I 


N ^ 1/2 

E / m tr(t^(RN)e(RN))^^x < / tr{cr{R^)e{RN))d\\ 

k=l'^^3£ ) \JClfj J 


<|6.28} then yields 


'Ll < Const ^ ^ tr(£r(RN)e(RN))dx'j ^ S(Yfc)dx^ 


VA :=1 •'% 


(6.29) 


6.2.3 Estimate for E 3 


Owing fo fhe Beffi formula. Lemma and fhe assumption fhaf fhe supporf of fhe cuf-off funcfion is 


confained in , we deduce 


r,_fr(n(RN)e(;cF{J(x - oW)i/;« - Y,}))dx 

•^®3£ \'^E 

If fhen follows fhaf 

^fr(tr(RN)e(;t«Y,))dx= ^tr((r(RN)e(;c^'^){Tfc + J(x - oW)i/;« - Y,}))dx . (6.30) 

The symmetry of the functional on the right-hand side implies 

/ ^fr(tr(RN)e(;cF{Y, - J(x - 0«)i/>« - Y^}))dx 

jBo \a;£ 


= ^fr(tr(xF{Y,+J(x-0«),/;«-Y,})e(RN))dx. 

J \(x)f 


(6.31) 


After applying the Cauchy and Schwarz inequalities to (6.311 and combining the result with ( 6.30) , it can be 
derived that 


y^(^^^tr(£r(4'^4Y,+J(x-OW)i/;W -Y,})n(4'^4Y,+J(x-OW)t/.« -Y,}))dxj (6.32) 


where S(U) is given in (2.6l. Then (6.17 1 and ( 6.18} provide 

j -^^ir{cr{RMxf'^^k))dx < Const (S(4'^4'r)t + J(x - oW)i/;W - Y^})d) 

\w^ \JB^^ 


1/2 


(6.33) 


Here, as a result of fhe inequalify 

S(wv) < Consf{|VM|^|v|^ + M^S(v)} 
for any vecfor funcfion v and scalar funcfion u, if can be asserfed fhaf 

+ J(x - 0«)i/>W - Y^})dx 

B'i^ 

f{e-2 |Yfc+J(x-OW)«/>« -Y^l^dxA f^^^S(^k)dx\ . 


°3e 

< Consf-j 
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(k) 

Again applying the Poincare inequality and the Friedrichs inequality in to the first integral on the above 
right-hand side (similar to <|6.24[ and 16.251) gives 


+ J(x - 

This estimate together with |6.33| yields 


X < Const 




^3 < Const 


/gW (/gW^^tr(£r(RN)e(RN))dx 


1/2 


6.2.4 Proof of ^6.11t 

Therefore 16.141, < 6.15^ , < 6.19} , ( 6.29} and (6.341 asserf fhaf 


N 


f 0,1 


fr(£r(Rf^)e(Riv))dx < Consf / S(;i(:oYo)dx + ^ S{Yi,)dx 


k=l •' % 


As a result of ( |2.6} , for a vector function v 

S(v) < Consf |E(v)|^ , 

and fhis wifh the definition of xo ( 6.35} yields ( 6.11} . 

6.3 Proof of Lemma and Theorem 1. Estimation of the energy for R 

The inequalify ( [6.11} leads to 

f tr(p'(Rf^f)e(RN))dx < Const {JC + C + M+J^} , 

J Dju 


where 


/C= / |Yo|2dx+ [ |E(Yo)|2dx, 

Jv Jv 

N . 

^ = E / w |E[u(x) - E(x - 0 «)(S(V,)^u(x)) 

k=l 

EL e( E {qF^(x)-(S(Vw)^H(w,x))^M, 

^ j^k 


x=0(9 


dx, 


M = 


(/) 


l<j<N 


v=o9') 


} c(^') 


- ^ S(x-OW)S(Vx)^(E(Vw)^G(w,x))^Mp')cW 


l<j<N 


x=0(9 

w= 0 // 


dx, 


.V=E/„5.,K(D(Vw)^H(w,x))WcI«|_^^,.,) 


dx . 


Owing to the representation of Yg in ([6.4} and Lemmathe term 1C admits the estimate 


1C < Const £' 


.8 1 


{ f (f 1A+ f (f )A} 

Wv ^|x - O0)|3/ Jv V ^ |x - 0(/)|4/ / 


|c0')| 


N 


/=1 


N 


N 


< Const 1^ Y] 1 /-rrz 

^tlUv|x-0(/)|6 


dx 

1 


|CW| 


Iv |x-O 0 )| 8 ^^} ' 


;= 

dx -|- 


(6.34) 


(6.35) 


(6.36) 

(6.37) 

(6.38) 


(6.39) 


(6.40) 
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where the last estimate has been obtained through the Cauchy inequality. Since dist(3n, 3<j;) > 1, the final 
estimate for 1C, affer applying Lemma is 

N 

1C < Consf e^d~^ < Consf . (6-41) 

M 

To esfimafe £, fhe Taylor approximafion is used fo expand fhe firsf-order derivatives of fhe function u abouf 
X = , as follows 



|E[u(x)]-E[u(x)] 




^dx 


N 

< Consf IIV (g) E[u(x))] 

P=i 



A local regularify esfimafe for fhe second-order derivafives of fhe componenfs of u inside co fhen leads fo 

£ < Consf e5d"3||E(u)|||^(^) . (6.42) 

By using fhe boundary condition for fhe regular parf H (see Section 3), fhe ferm A4 can be wriffen in fhe 
form 


N 

-M = E 


d(^) 




E [E(QW(x)-(E(Vw)^r(w,x))' 

j^k 


l<j<N 

-E ((S(Vw)^G(w,x))^ 






X 

II 

O 



v=o(i'i 


dx . 


Nexf, using Lemma and fhe Taylor expansion abouf x = of fhe second-order derivafives of fhe com¬ 
ponenfs of G, esfablishes fhe esfimafe 


N 


M < Consf E E 


|C0')| 


dx 


k=i j^k |OW-O0-)|4 

1<;<N 

N N 

< Consf e“ y |C(P)|2 T E ^ 

h k'k |oW-o(/)|8 

l<j<N 


1 


(6.43) 


Lemma 1^ fhen yields fhe final esfimafe for A4: 


N 


M < Consf E 


dxdy 

vxw: Iy _ ,,|S 

X—y|>rf I 


N 


< Consf E ^ Consf ||E(u) 


P=i 


\Loo(n] 


(6.44) 


Since fhe derivafives of fhe componenfs of H are bounded wifhin fhe cloud co, we deduce 


N 

J\f < Consf E |C^^^|^ < Consf e^d“^||E(u) ||E(n) • (6.45) 

k=l 

The energy esfimafe confained in Lemma|4|is fhen proved by combining ||6.41|, (6.42l, (6.44), (6.45) and (|6.36 1. 

□ 

Now we prove Theorem If remains fo consider fhe formal approximafion for U]v in Lemma which 
relies on fhe solvabilify of a parficular algebraic sysfem (0. The solvabilify of fhis sysfem was proved in 
Lemma which fogefher wifh fhe energy esfimafe in Lemma proves Theorem □ 


6.421, (6^, (6^ and (^ 
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7 Illustration: Simplified asymptotic formulae 

In this section, we present simplified asymptotic formulae for un in fhe far-field region away from fhe cloud 
of voids and also in fhe case when an infrnife elasfic medium confaining fhe cloud is considered. If is also 
shown in Appendix C fhaf for spherical voids, fhe model boundary layers of Problem 3 of Secfion can be 
consfrucfed explicifly in fhe closed form, along wifh fhe dipole mafrices for fhese spherical cavifies. 


7.1 Far-field approximation to Un 

Given fhe dipole mafrices 1 < k < N, fhe asympfofic formula {1.6 > of Theorem]^ is simplified under 
fhe consfrainf fhaf fhe poinf of measuremenf of fhe displacemenf is disfanf from fhe cloud of voids. 


Corollary 1 Let disf(x,a;) > 1. The asymptotic formula for un admits the form 


N 

un(x) = u(x) + ^(E(Vz)’^G(z,x))^M, 
k=l 


{k) 


s=oW 


cW+Fn(x), 


(7.1) 


where k = 1,...,N, satisfy the system (1.7), 


N 


Fn(x) = O ^ 


|cW| 


\k=l 


|x-OW|4 


R 


N / 


□ 


and satisfies (1.8). 

Proof. Formula < |7.1| follows from Lemma 

If is nofed fhaf in fhe simplified represenfafion < |7.1| for u^, informafion abouf fhe small voids is confained 
in fheir dipole characferisfics represenfed by , 1 < A: < N. In particular, if fhe voids are spherical cavifies 


of radius flj wifh cenfre 1 < k < N, fhen fhe dipole mafrix is given by 


Ik) ^ (A + 2p)7r(flf^)3 
y{9\ + Uy) 


03x3 

03x3 


(7.2) 


wifh 


= 


m m — 40//^ m — AQyf 
m — 
m 


m 


m — 40y^ 

m — 40y^ m — 40y^ 

= 40y%^3 . 


m = 9A^ + 20Ay + 36y^ 


(7.3) 


If is nofed fhaf fhe mafrix ' for fhe spherical cavify in fhe infinife space is negafive definife. Thus (7.21, 
< |7.3[ ), fogefher wifh Corollary]^ gives fhe far-field approximation for u^ in an elasfic solid confaining a cloud 
of arbifrary spherical cavifies. 


7.2 Far-field approximation for un in an infinite elastic medium with a cloud of voids 

Here we consider fhe problem when O = ]R^, so fhaf 0^ = is fhe infinife space confaining a 

cloud of voids. 

In fhis scenario, we search for fhe approximation fo un which is now a solution of fhe problem: 

h(Vx)uN(x)-Ff(x) = O , xen^, (7.4) 

T„(Vx)un(x) = O , X e dcoP , 1 < ; < N , 
un(x) = 0(|x|“^) , for |x| ^ 00 . 
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(7.5) 

(7.6) 














The vector function f is also supplied with the conditions that 




X X f(x)dx = O , 


and the support of f, as before, is chosen fo safisfy disf(da;,supp f) = 0(1). 

Finally, before sfafing resulfs concerning fhe approximafion of u^, we furfher infroduce some model 
quanfifies. We require fhe field u which solves fhe problem E) and < |3.2| , and fhaf is also supplied wifh fhe 
additional condition of decay af infinify <|7.6|. The mafrix 


P = 


S(Vx)^(E(Vy)^r(y,x))^ 
06x6 / 


x= 0 (') ' 
y= 0 (i) 


if > 

ofherwise, 


is also needed in fhe nexf resulf. We nofe fhaf in fhe considered case fhe regular parf H = 0, so fhaf Green's 
fensor in O is G = T, fhe Kelvin-Somigliana fensor, which is defined in ( |3.6| . 

Firsf, as a direcf consequence of Corollaryj^we have 

Corollary 2 Let dist{x,co) > 1, then the asymptotic formula for ujv admits the form 


N 


un{x) = u(x) + ^(S(Vz)^F(z,x))^M, 
k=l 




5 = 0 ('') 


CW +91n(x) 


(7.7) 


where 


(ff ehcWl 
mM = o[Y. 


vSlx-oW|4; 

C = ((Ci^i)^,..., (Ci^i)^)^ solves the linear algebraic system 

- V = (]l6Nx6N + PM)C , 


R 


N / 


(7.8) 


and Rjv/ satisfies (1.8). 


Once again, fhe dipole mafrix for a spherical cavify (see < |7.2| , | |7.3| ) can be used wifh 17.71 fo describe fhe 
far-field behaviour of u^ in an infrnife elastic space confaining a cloud of spherical cavities. 


7.3 Uniform approximation for in the infinite elastic space containing a cloud of 
voids 

Corollary 2 can be exfended fo a uniform approximafion xi]sj, safisfying 


7.41, (7.6 1 , inside 


O'). 


Corollary 3 Let the small parameters e and d satisfy the inequality 


e < cd, 

where c is a sufficiently small constant. Then the approximation for u^ is given by 

Un(x) = u(x) + qF(x)c(^) + Rn(x) , 

J:=l 


(7.9) 


(7.10) 


and Rj>j satisfies (1.8). 


ik) 

Mafrices such as Q) ' can be consfrucfed in fhe explicif closed form for cerfain geomefries. For spherical 
voids, fhe represenfafion of fhis mafrix is given in Appendix C. Thus, if fhe cloud co is composed of a non¬ 
periodic arrangemenf of spherical voids 1 < j < N, fhen fhe approximafion sfafed in fhe previous 
Corollary, fogefher wifh fhe represenfafion of fhe mafrix in Appendix C is readily applicable here. 
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8 Concluding remarks 


A uniform asymptotic representation for a solution of a mixed boundary value problem of elasticity has been 
constructed and justified for a solid containing a cloud of many voids. This extends significantly the results 
of the papers IIT8ll22ll24l on meso-scale asymptotic approximations of fields in domains with multiple defects. 
It is worth noting that the asymptotic representation < |1.6| of Theorem 1, contains an important information 
about the dipole fields of a meso-scale cloud of voids. In addition to the sum of individual contributions from 
the dipole fields of small voids, we have also obtained a term characterising a mutual interaction between 
the voids, which is often neglected in the dilute approximation procedures. This result is significant in the 
area of applications linked to non-destructive testing of porous solids, where a position of a cloud and its 
composition can be identified through the use of the asymptotic formula ( 1.61 accompanied by the boundary 
measurements for different test loading conditions. 

It is also essential to note that the meso-scale approximation | |1.6[ l is valid for different shapes of small 
voids when e < Const d, for a sufficiently small constant, and this surpasses the range of applicability of the 
homogenisation approximations. 
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Appendix A: Local regularity of solutions to the homogeneous Lame sys¬ 
tem 


Here, a result concerning the estimate for the derivatives of fhe solufion fo fhe homogeneous Lame sysfem 
via their anti-derivatives is derived: 

Lemma 5 Let w = be a solution of the homogeneous Lame system in a domain O and let Br c O, with 

Br = {x: |x| < _R}, then the estimate 


^^(O) < Const R ^suplwl 


(A.l) 


holds. 


The proof of fhe lasf esfimafe uses fhe mean value fheorem for vecfor funcfions safisfying fhe homoge¬ 
neous Lame sysfem, as discussed below and in Il4l . 

Lemma 6 Let w = be a solution of the homogeneous Lame system in a domain O and Br c O, with 

Br = {x : |x| < R}, then: 


(i) 


(ii) 


U - 8nRHA + {p) U 


3(A-f) 


Wi(O) = 


75{X + y) 
87 rR 5 (A-|- 4 /r) Jbr 


/ XiXjWj{x)dx — 
^ Bj? 


15(A-f) f 

8nR^{A + Ay) Jbr 


^Wi{x)dx . 


(A.3) 


Proof, (i) The mean value fheorem of I A.2| was proved in |(4|. (ii) To derive < A.3) , apply < A.2) inside fhe 
ball Br C n. Then mulfiplying fhrough fhe resulting equation by and integrating both sides with respect 
to r between zero and R yields <|A.3 1. 


Proof of (A.l): The mean value theorem I A.3 1 is applied in Br to the function as follows: 

aXk 


dwi 


75{\ + y) 

dxk 8nR^{X +Ay) Jbr ' ^ dx^^ 


x)dx — 


15{A-y) I ^ _ 
87rR5(A-|-4;i) Jbr dx/^ 




(x)dx . 


Infegration by parfs fhen yields fhe fwo identities: 

dw 


XiXj^—^(x)dx = — 

idx^ 


j (SikXjWj + XiW}f)dx -|- J n^XiXjWjdsx , 


f 7 f 


Then <|A.5) and ||A. 6 ) give the estimates: 


ibr 


r dw; , 

/ XiXj^^{x)dx < Const R sup |w| and 
Ibr ’ dXk 


ni^\x\^Wjdsx . 


IBBr 

[ IxP—^(x)dx < Const R^ sup |w| 

IBr dXfc Bk 


and combining these with |A.4l yields the local regularity estimate <|A.l). The proof is complefe. 


□ 


(A.4) 

(A.5) 

(A. 6 ) 

(A.7) 

□ 


Appendix B: Proof of (|5.6|) 


Here, fhe proof of (|5.6) is carried ouf by firsf developing an identify which will lead fo an infegral represen- 
fafion of ||5.5) in Section B.l Then we prove some auxiliary infegral identifies in Section B.2 fhaf are used fo 


complefe fhe proof of <|5.6) in Section B.3 
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B.l Poisson-type representation of the second-order derivatives of Green's tensor 

The proof of fhe nexf lemma uses fhe mean value fheorem for solutions of fhe homogeneous Lame sysfem 
inside disjoinf balls denofed by = {x : |x — \ < d/4}, j = 1,... ,N. 

According fo ID and Lemma 10 of Il24ll , fhe nexf resulf holds: 


Lemma 7 For j ^ k, 1 < i,k < N, the identity 

Ssp (Vz) E/jg (Vw) Gs/, (Z, W) 

is valid, where 


w 


(87r)2(f)6(A + 4^)2 


:A 


(/A) 


(B.l) 


Ap'q ^ — 36^(A + 4p)^/7’p^''^' ^ + 90(A+ 4/r)(A + ^ 

- 18(A + 4;^)(A - y) + 225(A + yfj!;, 


-45(A^ - y^) [ + 9(A - 


'v<t 


(B.2) 

(B.3) 


and the terms for 1 < s < 6are 


Jpq — 

Li) . 


^(2//A) _ 

Jpc; — 

Li) . 

Lk) 

Jpq — 

Li) . 

Lk) 

7" (4/A) _ 

Jpq — 

Li) . 

fsik) 

^(5,j,k) _ 
Jpq — 

Li). 

Lk) 

Jpq — 

Li) ■ 

fsik) 


JbH) /bW S«p(^z)SN(^w)G,fc(Z, W) dWdZ , 

Li) /bW ' 

- oW),^(E,p(Vz)Ej,g(Vw)G„b(Z, W)) dWdZ , 

- oW)t(W- 


.gZ^3^j,vS«p(Vz)Efe,(Vw)Gts(Z,W))dWdZ, 
,g^^(E„,p(Vz)E„q(Vw)Gfb(Z,W))dWdZ, 


L) L)~ ~ o‘'^)^a^;^(s.p(Vz)s,,(Vw)G„,(z, w)) dwdz. 

Before presenting fhe proof, if is nofed fhaf | |B.1| is also a cormecfed wifh fhe classical resulfs of ||TJ |2l on 
esfimafes for solutions of ellipfic partial differential equations (eg. see Theorems 7.3 in |[T| and Theorem 9.3 
in lEI). 

Proof. Using the Kronecker delfa fhe ferm in fhe leff-hand side of is rewritten as 


^ms^hn'^mp (Vz) '^nq {Vw)Gsh{Z,W) 


z=o)i) ■ 
W=o(9 


(B.4) 


From here, fhe ferm Emp(Vz)S„^(Vw)Gs;j(Z, W) 2 -o(;) considered as enfries of fhe mafrix 

W=0(9 

Smp(Vz)E„^(Vw)G(Z, W) 

which safisfies fhe homogeneous Lame equafion for Z G B^if The franspose of fhe preceding mafrix also 
satisfies fhe homogenous Lame equafion for W G B^^fkf^j. Then rep eatin g the steps of fhe proof of Le mma 
10 in 1241 , for fhe lasf mafrix and inserting fhe resulfrng expression in I B.4 1 we arrive af Apq^'^ (see (B.3 ) and 
fhe relation B- The proof is complefe. □ 
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B.2 Auxiliary integral identities 

Now the Poisson-type representations for the second-order derivatives of Green's fensor are in place, furfher 
idenfifies are now derived which are used in fhe proof Lemma in fhe next section. From here, we will also 
make use of fhe vecfor and mafrix funcfions 


and 


respecfively. 

Lemma 8 The identities 


hold. 


«I»(x) 


Mp'^cW, ifxesW, 
O , ofherwise, 


0(x) 


® (x - OW) , if x e BW , 
O, ofherwise. 


P 3 

/ 0,„t(Z) —(E„p(Vz)Gfb(Z,W))dZ = 0, 

Jo, dZjfi 

cl>p(z)E«p(Z)Gi„(W,Z) dZ = 0, 


(B.5) 


(B.6) 


(B.7) 

(B.8) 


Proof. We prove izi and nofe fhaf fhe idenfify | |B.8| is proved in a similar way wifh obvious modificafions. 
Sef 

/(W)= / J^E,p(Vz)G(W,Z)0^(Z)dZ, 

Jn oZ.„ 

which is fhe same as fhe leff-hand side in < |B.7| , wifh fhe assumpfion fhaf fhe subscripf indices in fhe above 
are free. The funcfion / is fhen a 3 x 6 mafrix whose columns satisfy fhe homogeneous Lame sysfem. Indeed, 
since affer an application of fhe Lame operator, it is possible to retrieve through the definition of G: 

-L(Vw)/(W) = / ^(E,p(Vz)<5(W-Z)l3)0^(Z)dZ 
Jci aZn 

= [ s{z-w)^Zap{Vz)e^{z)dz. 

Jci aZn 

Now, when considering fhe cases W G and W G n\ BO), the definition of 0 shows fhat fhe 

above righf-hand side is equal fo 03 x 6 - 

Again the definition of G also ensures fhat /(W) = © 3 x 6 for W G 90. An application of fhe Betti's 
formula fo /(W) and Green's mafrix G in O shows fhen shows fhaf /(W) = © 3 x 6 for W G O and fhe proof 
of <|B.7| is complefe. □ 


B.3 The estimate for ( |5.5^ 


Relation (5.6 1 is then a result of the next Lemma: 


Lemma 9 The relation 


is valid, where 


(MC,SMC) = - 


(87r)2(|)6(A + 4^0 


2 ' 


N 


/=l 

with repeated subscript indices being regarded as the indices of summation and 

|(MC,SMC)| < Gonstd"^(MC,MC) . 


(B.9) 


(B.IO) 


(B.ll) 
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Proof. Representations (B.9) and (B.IO). The combination of IB.I I and 1 5.5 1 then delivers the expression 


(MC,SMC) = 


g-h 


i8nm)HA + Ap) 


2 ' 


with 


N N 

EE 

y=l k=l 




and h is defined in I B.IO^ . 

The above expression for g admits a form that makes use of (|B.5| and ||B.6| 


g = + +180(A + 4fi)(A + p)/C(2) -36(A + 4/r)(A-fi)/C(3) 

+225(A + ?i)2/cW - 90 (a2 - + 9(A - , 


(B.12) 


(B.13) 


(B.14) 


where 

/ / cl>p(z)cI),(W)E„p(Vz)Efc<,(Vw)G„j,(Z,W)dWdZ, 

/c(2)= [ / 0pt(Z)cI)^(W)^(E«p(Vz)Efc<,(Vw)Gtfe(Z,W))dWdZ, 

JCi JCl 

/C(3)=/ [ 0p,(Z)cI>^(W)^(E«p(Vz)Efcq(Vw)G„b(Z,W))dWdZ, 

JCi JCl C/Zis 

= J^J^&pt{Z)Qc!s{W)^^^^{Eap{Vz)Eh^{Vv,)Gts{Z,W))dWdZ, 
= |^|^0p«(Z)0<,,(W)^^^(E„p(Vz)Efc,(Vw)G„j,(Z,W))dWdZ. 


Now, the term is rewritten using the Kronecker delta as 

/c(2) = 


In 


f 3 

/n O^n 


'ap 


(Vz)Gtfc(Z,W))dZdW, 


where as shown in Lemmathe irmer integral is zero. Thus = 0. Similar conversions and Lemma 
also show that the terms K.^, lC^i\ 3 < j < 6 are equal to zero. In this way, we have shown that ^ = 0 (see 
< B.13| ). The proof of |B.9l and | B.10| i s complete. 

Estimate for h. Next, to prove (5.6 1 , the estimate for the quantity h 

\h\ < Const d3(MC,MC) , 


is proved. 

To show this, an estimate for the terms 1 < m < 6 is needed that make use of the fact that for 

x,y e n, ||G(x,y)|| = 0(|x — y|“^). Employing this, a majorant for is given by 

< Const < Const d=>. (B.15) 

The estimates for 2 < m < 6 are developed in a similar way. Next, consider the term 

N 

,1 < m < 6 . 

/=i 

Recalling that subscript indices are the indices of summation, and repeatedly applying the Cauchy inequality 
the above admits the inequality 

E(mF>cW),j,',"«'(mS''c(»), < E imF’cwp ( E 

7=1 7=1 \m=1 
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for 1 < m < 6. The preceding combines with lB.15 1 and the estimates for 2 < m < 6 to show that 


N . ... N . 

< Const 

M i=i 


Therefore, consulting | |B.9| it can be asserted that | B.11| holds. Thus the proof of the present lemma and < |5.6| 
is complete. □ 


Appendix C: Explicit representation of dipole fields for spherical cavities 

It is shown in this section that for certain geometries, model fields used in the asymptotic approximations 
presented here can be constructed in the closed form. Here, it is assumed that the voids co^J\ j = 1,... ,N, are 
spherical cavities. The matrix for a spherical cavity with radius and centre at = {O® 
in an infinite solid can be reconstructed using the approach presented in IIT^ that makes use of Papkovich- 
Neuber potential representation for solutions to three-dimensional elasticity problems. 

In this case, the matrix takes the form 

qF(x) = -(E(Vx)^r(x,oW))^MF + ^_l^E(x - 1^{e(x - o«)2i(x - o«) 

-h(xi -0®)(X2 -Of^)(x3 -hS[ll(x-OW)E(x-OW)2t®} . 

(C.l) 


Here the dipole matrix is given in (7.21, (7.31, and the matrices 1 < p < 3 are 


(;t)_ 3(A + p)(flf)5 


21} ' = - 


21} ^ = 


9A 14p 


9A 14p 


03 ' 

, 23 ( 1 ) = 

■3 1 r 

13 1 

_ U 3 2 JI 3 


1 1 3_ 

- [O 3 25 ( 2 ) 

, = 

'0 0 1‘ 
0 1 0 



1 0 0 


(k) _ 30(A + p)(af))5 
9A 14p 


2ir = 


03 03 

03 13 


Also the matrix functions in (C.l I, 2) and SJl, are given as 


m{x) = 


and 


21(x) = 


15(A + p)(fl®)5 
9A 14p 


2 )(x) O 3 
O 3 O 3 


xf 0 0 

0 x} 0 

0 0 x^j 


®(x) = 


'xf 

x} 

4' 


'1 

1 

1' 

x} 

x} 

4 

= 

1 

1 

1 

x^ 

x} 

4\ 


1 

1 

1 


S[)l(x) . 
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